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METHODS OF PREDICTING MUSCULOSKELETAL DISEASE 
Technical Field 

[0001] This invention relates to using imaging methods for diagnosis, 
prognostication, monitoring and management of disease, particularly where 
5 that disease affects the musculoskeletal system. This invention identifies 
novel imaging markers for use in diagnosis, prognostication, monitoring and 
management of disease, including musculoskeletal disease. 

Background 

[0002] Osteoporosis and osteoarthritis are among the most common 
10 conditions to affect the musculoskeletal system, as well as frequent causes of 
locomotor pain and disability. Osteoporosis can occur in both human and 
animal subjects (e.g. horses). Osteoporosis (OP) and osteoarthritis (OA) 
occur in a substantial portion of the human population over the age of fifty. 
The National Osteoporosis Foundation estimates that as many as 44 million 
15 Americans are affected by osteoporosis and low bone mass. In 1997 the 
estimated cost for osteoporosis related fractures was $13 billion. That figure 
increased to $17 billion in 2002 and is projected to increase to $210-240 
billion by 2040. Currently it is expected that one in two women over the age 
of 50 will suffer an osteoporosis-related fracture. 

20 [0003] Imaging techniques are important diagnostic tools, particularly for 
bone related conditions such as OP and OA. Currently available techniques 
for the noninvasive assessment of the skeleton for the diagnosis of 
osteoporosis or the evaluation of an increased risk of fracture include dual x- 
ray absorptiometry (DXA) (E.astell et al. (1998) New Engl J. Med 338:736- 

25 746); quantitative computed tomography (QCT) (Cann (1988) Radiology 
166:509-522); peripheral DXA (pDXA) (Patel et al. (1999) J Clin Densitom 
2:397-401); peripheral QCT (pQCT) (Glueret. al. (1997) Semin Nuc! Med 
27:229-247); x-ray image absorptiometry (RA) (Glueret. al. (1997) Semin 
Nucl Med 27:229-247); and quantitative ultrasound (QUS) (Njeh et al. 

30 "Quantitative Ultrasound: Assessment of Osteoporosis and Bone Status" 
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1999, Martin-Dunitz, London England; U.S. Patent No. 6,077,224, 
incorporated herein by reference in its entirety). (See, also, WO 9945845; 
WO 99/08597; and U.S. Patent No. 6,246,745). 

[0004J DXA of the spine and hip has established itself as the most widely 
5 used method of measuring BMD. Tothiil, P. and D.W. Pye, (1992) Br J Radiol 
65:807-81 3. The fundamental principle behind DXA is the measurement of 
the transmission through the body of x-rays of 2 different photon energy 
levels. Because of the dependence of the attenuation coefficient on the 
atomic number and photon energy, measurement of the transmission factors 

10 at 2 energy levels enables the area densities (i.e., the mass per unit projected 
area) of 2 different types of tissue to be inferred. In DXA scans, these are 
taken to be bone mineral (hydroxyapatite) and soft tissue, respectively. 
However, it is widely recognized that the accuracy of DXA scans is limited by 
the variable composition of soft tissue. Because of its higher hydrogen 

1 5 content, the attenuation coefficient of fat is different from that of lean tissue. 
Differences in the soft tissue composition in the path of the x-ray beam 
through bone compared with the adjacent soft tissue reference area cause 
errors in the BMD measurements, according to the results of several studies. 
Tothiil, P. and D.W. Pye, (1992) Br J Radiol, 65:807-813; Svendsen, O.L., et 

20 al., (1995) J Bone Min Res 10:868-873. Moreover, DXA systems are large 
and expensive, ranging in price between $75,000 and $150,000. 

[0005] Quantitative computed tomography (QCT) is usually applied to 
measure the trabecular bone in the vertebral bodies. Cann (1988) Radiology 
1 66:509-522. QCT studies are generally performed using a single kV setting 

25 (single-energy QCT), when the principal source of error is the variable 
composition of the bone marrow. However, a dual-kV scan (dual-energy 
QCT) is also possible. This reduces the accuracy errors but at the price of 
poorer precision and higher radiation dose. Like DXA, however, QCT are 
very expensive and the use of such equipment is currently limited to few 

30 research centers. 
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[0006] Quantitative ultrasound (QUS) is a technique for measuring the 
peripheral skeleton. Njeh et at. (1997) Osteoporosis Int 7:7-22; Njeh et at. 
Quantitative Ultrasound: Assessment of Osteoporosis and Bone Status. 1999, 
London, England: Martin Dunitz. There is a wide variety of equipment 
5 available, with most devices using the heel as the measurement site. A 
sonographic pulse passing through bone is strongly attenuated as the signal 
is scattered and absorbed by trabecular Attenuation increases linearly with 
frequency, and the slope of the relationship is referred to as broadband 
ultrasonic attenuation (BUA; units: dB/MHz). BUA is reduced in patients with 

10 osteoporosis because there are fewer trabeculae in the calcaneus to 

attenuate the signal. In addition to BUA, most QUS systems also measure 
the speed of sound (SOS) in the heel by dividing the distance between.the 
sonographic transducers by the propagation time (units: m/s). SOS values 
are reduced in patients with osteoporosis because with the loss of mineralized 

15 bone, the elastic modulus of the bone is decreased. There remain, however, 
several limitations to QUS measurements. The success of QUS in predicting 
fracture risk in younger patients remains uncertain. Another difficulty with 
QUS measurements is that they are not readily encompassed within the WHO 
definitions of osteoporosis and osteopenia. Moreover, no intervention 

20 thresholds have been developed. Thus, measurements cannot be used for 
therapeutic decision-making. 

[0007] There are also several technical limitations to QUS. Many devices 
use a foot support that positions the patient's heel between fixed transducers. 
Thus, the measurement site is not readily adapted to different sizes and 
25 shapes of the calcaneus, and the exact anatomic site of the measurement 
varies from patient to patient. It is generally agreed that the relatively poor 
precision of QUS measurements makes most devices unsuitable for 
monitoring patients' response to treatment. Gluer (1997) J Bone Min Res 
12:1280-1288. 

30 [0008] Radiographic absorptiometry (RA) is a technique that was 

developed many years ago for assessing bone density in the hand, but the 
technique has recently attracted renewed interest. Gluer et at. (1997) Semin 
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Nucl Med 27:229-247. With this technique, BMD is measured in the 
phalanges. The principal disadvantage of RA of the hand is the relative lack 
of high turnover trabecular bone. For this reason, RA of the hand has limited 
sensitivity in detecting osteoporosis and is not very useful for monitoring 
5 therapy-induced changes. 

[0009] Peripheral x-ray absorptiometry methods such as those described 
above are substantially cheaper than DXA and QCT with system prices 
ranging between $15,000 and $35,000. However, epidemiologic studies have 
shown that the discriminatory ability of peripheral BMD measurements to 

10 predict spine and hip fractures is lower than when spine and hip BMD 

measurements are used. Cummings et al. (1993) Lancet 341:72-75; Marshall 
et al. (1996) Br Med J 312:1254-1259. The main reason for this is the lack of 
trabecular bone at the measurement sites used with these techniques. In 
addition, changes in forearm or hand BMD in response to hormone 

15 replacement therapy, bisphosphonates, and selective estrogen receptor 

modulators are relatively small, making such measurements less suitable than 
measurements of principally trabecular bone for monitoring response to 
treatment. Faulkner (1 998) J Clin Densitom 1 -.279-285; Hoskings et al. 
(1998) N Engl J Med 338:485-492. Although attempts to obtain information 

20 on bone mineral density from dental x-rays have been attempted (See, e.g., 
Shrout et al. (2000) J. Periodonol. 71:335-340; Verhoeven et al. (1998) Clin 
Oral Implants Res 9(5):333-342), these have not provided accurate and 
reliable results. 

[0010] Furthermore, current methods and devices do not generally take 
25 into account bone structure analyses. See, e.g., Ruttimann et al. (1992) Oral 
Surg Oral Med Oral Pathol 74:98-1 10; Southard & Southard (1992) Oral Surg 
Oral Med Oral Pathol 73:751-9; White & Rudolph, (1 999) Oral Surg Oral Med 
Oral Pathol Oral Radiol Endod 88:628-35. 

[0011] The present invention discloses novel methods and techniques for 
30 predicting musculoskeletal disease, particularly methods and compositions 
that result in the ability to obtain accurate predictions about disease based on 
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bone mineral density and/or bone structure information obtained from images 
(e.g., radiographic images) and data. 

Summary of the Invention 

[0012] The invention discloses a method for analyzing at least one of bone 
5 mineral density, bone structure and surrounding tissue. The method typically 
comprises: (a) obtaining an image of a subject; (b) locating a region of interest 
on the image; (c) obtaining data from the region of interest; and (d) deriving 
data selected from the group of qualitative and quantitative from the image 
data obtained at step c. 

10 [0013] A system is also provided for predicting a disease. Any of these 
systems can include the steps of: (a) obtaining image data of a subject; (b) 
obtaining data from the image data wherein the data obtained is at least one 
of quantitative and qualitative data; and (c) comparing the at least one of 
quantitative and qualitative data in step b to at least one of: a database of at 

15 [east one of quantitative and qualitative data obtained from a group of 
subjects; at least one of quantitative and qualitative data obtained from the 
subject; and at least one of a quantitative and qualitative data obtained from 
the subject at time Tn. 

[0014] In certain aspects, described herein are methods of diagnosing, 
20 monitoring and/or predicting bone or articular disease (e.g., the risk of 

fracture) in a subject, the method comprising the steps of: determining one or 
more micro-structural parameters, one or more macroanatomicai parameters 
or biomechanical parameters of a joint in said subject; and combining at least 
two of said parameters to predict the risk of bone or articular disease. The 
25 micro-structural, macroanatomicai and/or biomechanical parameters may be, 
for example, one or more of the measurements/parameters shown in Tables 
1, 2 and/or 3. In certain embodiments, one or more micro-structural 
parameters and one or more macro-anatomical parameters are combined. In 
other embodiments, one or more micro-structural parameters and one or 
30 more biomechanical parameters are combined. In further embodiments, one 
or more macroanatomicai parameters and one or more biomechanical 
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parameters are combined. In still further embodiments, one or more 
macroanatomicaf parameters, one or more micro-structural parameters and 
one or more biomechanics parameters are combined. 

[0015] In any of the methods described herein, the comparing may be 
5 comprise univariate, bivariate and/or multivariate statistical analysis of one or 
more of the parameters. In certain embodiments, the methods may further 
comprise comparing said parameters to data derived from a reference 
database of known disease parameters. 

[0016] In any of the methods described herein, the parameters are 
10 determined from an image obtained from the subject, in certain 

embodiments, the image comprises one or more regions of bone (e.g., 
patella, femur, tibia, fibula, pelvis, spine, etc). The image may be 
automatically or manually divided into two or more regions of interest. 
Furthermore, in any of the methods described herein, the image may be, for 
15 example, an x-ray image, a CT scan, an MRI or the like and optionally 
includes one or more calibration phantoms. 

[0017] In any of the methods described herein, the predicting includes 
performing univariate, bivariate or multivariate statistical analysis of the 
analyzed data and referencing the statistical analysis values to a fracture risk 
20 model. Fracture risk models can comprise, for example, data derived from a 
reference database of known fracture loads with their corresponding values of 
macro-anatomical, micro-anatomical parameters, and/or clinical risk factors. 

[0018] In another aspect, the invention includes a method of determining 
the effect of a candidate agent on a subject's prognosis for musculoskeletal 

25 disease comprising: predicting a first risk of musculoskeletal disease in 
subject according to any of the predictive methods described herein; 
administering a candidate agent to the subject; predicting a second risk of the 
musculoskeletal disease in the subject according to any of the predictive 
methods described herein; and comparing the first and second risks, thereby 

30 determining the effect of the candidate on the subject's prognosis for the 
disease. In any of these methods, the candidate agent can be administered 
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to the subject in any modality, for example, by injection (intramuscular, 
subcutaneous, intravenous), by oral administration {e.g., ingestion), topical 
administration, mucosal administration or the like. Furthermore, the candidate 
agent may be a small molecule, a pharmaceutical, a biopharmaceutical, an 
5 agropharmaceuticals and/or combinations thereof. 

£0019] In other aspects, the invention includes a kit that is provided for 
aiding in the prediction of musculoskeletal disease (e.g., fracture risk). The kit 
typically comprises a software program that uses information obtained from 
an image to predict the risk or disease (e.g., fracture). The kit can also 
10 include a database of measurements for comparison purposes. Additionally, 
the kit can include a subset of a database of measurements for comparisons. 

[0020] In any of these methods, systems or kits, additional steps can be 
provided. Such additional steps include, for example, enhancing image data. 

[0021] Suitable subjects for these steps include for example mammals, 
15 humans and horses. Suitable anatomical regions of subjects include, for 
example, dental, spine, hip, knee and bone core x-rays. 

[0022] A variety of systems can be employed to practice the inventions. 
Typically at least one of the steps of any of the methods is performed on a 
first computer. Although, it is possible to have an arrangement where at least 
20 one of the steps of the method is performed on a first computer and at least 
one of the steps of the method is performed on a second computer. In this 
scenario the first computer and the second computer are typically connected. 
Suitable connections include, for example, a peer to peer network, direct link, 
intranet, and internet. 

25 [0023] It is important to note that any or all of the steps of the inventions 
disclosed can be repeated one or more times in series or in parallel with or 
without the repetition of other steps in the various methods. This includes, for 
example repeating the step of. locating a region of interest, or obtaining image 
data. 
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[0024] Data can also be converted from 2D to 3D to 4D and back; or from 
2D to 4D. Data conversion can occur at multiple points of processing the 
information. For example, data conversion can occur before or after pattern 
evaluation and/or analysis. 

5 [0025] Any data obtained, extracted or generated under any of the 
methods can be compared to a database, a subset of a database, or data 
previously obtained, extracted or generated from the subject. For example, 
known fracture load can be determined for a variety of subjects and some or 
all of this database can be used to predict fracture risk by correlating one or 
10 more macro-anatomical or structural parameters (Tables 1, 2 and/or 3) with 
data from a reference database of fracture load for age, sex, race, height and 
weight matched individuals. 

[0026] The present invention provides methods that allow for the analysis 
of bone mineral density, bone and/or cartilage structure and morphology 

15 and/or surrounding tissue from images including electronic images and, 
accordingly, allows for the evaluation of the effect(s) of an agent (or agents) 
on bone and/or cartilage. It is important to note that an effect on bone and/or 
cartilage can occur in agents intended to have an effect, such as a therapeutic 
effect, on bone and/or cartilage as well as agents intended to primarily effect 

20 other tissues in the body but which have a secondary, or tangential , effect on 
bone and/or cartilage. The images (e.g., x-ray images) can be, for example, 
dental, hip, spine or other radiographs and can be taken from any mammal. 
The images can be in electronic format. 

[0027] The invention includes a method to derive quantitative information 
25 on bone structure and/or bone mineral density from an image comprising (a) 
obtaining an image, wherein the image optionally includes an external 
standard for determining bone density and/or structure; and (b) analyzing the 
image obtained in step (a) to derive quantitative information on bone 
structure. The image is taken of a region of interest (ROI). Suitable ROI 
30 include, for example, a hip radiograph or a dental x-ray obtained on dental x- 
ray film, including the mandible, maxilla or one or more teeth. In certain 

8 



WO 2004/062495 



PCT/US2004/000315 



embodiments, the image is obtained digitally, for example using a selenium 
detector system, a silicon detector system or a computed radiography system. 
In other embodiments, the image can be digitized from film, or another 
suitable source, for analysis. 

5 [0028] A method is included where one or more candidate agents can be 
tested for its effects on bone. Again, the effect can be a primary effect or a 
secondary effect. For example, images obtained from the subject can be 
evaluated prior to administration of a candidate agent to predict the risk of 
disease in the absence of the agent. After administration of the candidate 

10 agent(s), an electronic image of the same portion of a bone of the subject can 
be obtained and analyzed as described herein to predict the risk of 
musculoskeletal disease. The risk of disease prior to administration of the 
candidate agent and after administration can then be compared to determine 
if the agent had any effect on disease prognosis. Information on bone 

15 structure can relate to a variety of parameters, including the parameters 
shown in Table 1 , Table 2 and Table 3, infra. The images or data may also 
be compared to a database of images or data (e.g., "known" images or data). 
The candidate agent can, for example, be molecules, proteins, peptides, 
naturally occurring substances, chemically synthesized substances, or 

20 combinations and cocktails thereof. Typically, an agent includes one or more 
drugs. Further, the agent can be evaluated for the ability to effect bone 
diseases such as the risk of bone fracture {e.g., osteoporotic fracture). 

[0029] In any of the methods described herein, the analysis can comprise 
using one or more computer programs (or units). Additionally, the analysis 

25 can comprise identifying one or more regions of interest (ROI) in the image, 
either prior to, concurrently or after analyzing the image, e.g. for information 
on bone mineral density and/or bone structure. The bone density information 
can be, for example, areas of highest, lowest or median density. Bone 
structural information can be, for example, one or more of the parameters 

30 shown in Table 1, Table 2 and Table 3. The various analyses can be 

performed concurrently or in series. Further, when using two or more indices 
each of the indices can be weighted equally or differently, or combinations 
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thereof where more than two indices are employed. Additionally, any of these 
methods can also include analyzing the image for bone mineral density 
information using any of the methods described herein. 

{0030] Any of the methods described herein can further comprise applying 
5 one or more correction factors to the data obtained from the image. For 
example, correction factors can be programmed into a computer unit. The 
computer unit can be the same one that performs the analysis of the image or 
can be a different unit. In certain embodiments, the correction factors account 
for the variation in soft-tissue thickness in individual subjects. 

10 [0031] These and other embodiments of the subject invention will readily 
occur to those of skill in the art in light of the disclosure herein. 

Brief Description of the Drawings 

[0032] Figs. 1a and b are block diagrams showing the steps for extracting 
data from an image and then deriving quantitative and/or qualitative data from 
15 the image. 

[0033] Figs. 2a-c are diagrams showing an image taken of a region of 
anatomical interest further illustrating possible locations of regions of interest 
for analysis. 

[0034] Figs. 3a-j illustrate various abnormalities that might occur including, 
20 for example, cartilage defects, bone marrow edema, subchondral sclerosis, 
osteophytes and cysts. 

[0035] Figs. 4a and b are block diagrams of the method of Fig. 1a showing 
that the steps can be repeated. 

[0036] Figs. 5a-e are block diagrams illustrating steps involved in 
25 evaluating patterns in an image of a region of interest. 

[0037] Fig. 6a-e are block diagrams illustrating steps involved in deriving 
quantitative and qualitative data from an image in conjunction with 
administering candidate molecules or drugs for evaluation. 
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[0038] Figs. 7a-d are block diagrams illustrating steps involved in 
comparing derived quantitative and qualitative information to a database or to 
information obtained at a previous time. 

[0039] Figs. 8a-d are block diagrams illustrating steps involved in 
5 comparing converting an image to a pattern of normal and diseased tissue 

[0040] Fig. 9 is a diagram showing the use one or more devices in the 
process of developing a degeneration pattern and using a database for 
degeneration patterns. 

[0041] Fig. 1 0 depicts regions of interest (ROIs) analyzed in Example 1 . 

10 [0042] Fig. 1 1 depicts results of biomechanical testing of 15 cadaveric hips 
and femurs. 

[0043] Fig. 12A-B, are reproductions of x-ray images depicting an 
exemplary induced fracture in cadaveric femur resulting from biomechanical 
testing and load. 

15 ' [0044] Fig. 13 is a graph depicting correlation of DXA femoral neck bone 
mineral density (BMD) versus biochemical fracture load as evaluated in 15 
fresh cadaveric hip samples. 

[0045] Fig. 14A-C are graphs depicting correlation of bone structure versus 
mechanical fracture load. Fig. 14A depicts correlation of maximum marrow 
20 spacing v. fracture load. Fig. 14B depicts correlation of maximum marrow 
spacing (log) v. fracture load. Fig. 14C depicts correlation of percentage of 
trabecular area v. fracture load. 

[0046] Fig. 15A-C are graphs depicting correlation of macro-anatomical 
features versus biomechanical fracture load. Fig. 15A depicts correlation of 
25 cortical thickness v. fracture load. Fig. 15B depicts correlation of hip axis 

length (HAL) v. fracture load. Fig. 15C depicts correlation of cortical thickness 
(standard deviation) versus fracture load. 
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[0047] Fig. 16 is a graph depicting multivariate analysis using a combination 
of bone structural and macro-anatomical parameters and shows the 
correlation of predicted fracture load to actual fracture load. 



5 Detailed Description 

[0048] The following description is presented to enable any person skilled 
in the art to make and use the invention, Various modifications to the 
embodiments described will be readily apparent to those skilled in the art, and 
the generic principles defined herein can be applied to other embodiments 

10 and applications without departing from the spirit and scope of the present 
invention as defined by the appended claims. Thus, the present invention is 
not intended to be limited to the embodiments shown, but is to be accorded 
the widest scope consistent with the principles and features disclosed herein. 
To the extent necessary to achieve a complete understanding of the invention 

15 disclosed, the specification and drawings of all issued patents, patent 

publications, and patent applications cited in this application are incorporated 
herein by reference. 

[0049] The practice of the present invention employs, unless otherwise 
indicated, currently conventional methods of imaging and image processing 

20 within the skill of the art. Such techniques are explained fully in the literature. 
See, e.g., WO 02/22014, X-Ray Structure Determination: A Practical Guide, . 
2 nd Edition, editors Stout and Jensen, 1989, John Wiley & Sons, publisher; 
Body CT: A Practical Approach, editor Slone, 1999, McGraw-Hill publisher; 
The Essential Physics of Medical Imaging, editors Bushberg, Seibert, 

25 Leidholdt Jr & Boone, 2002, Lippincott, Williams & Wilkins; X-ray Diagnosis: A 
Physician's Approach, editor Lam, 1998 Springer-Verlag, publisher; Dental 
Radiology: Understanding the X-Ray Image, editor Laeb'tia Brocklebank 1997, 
Oxford University Press publisher; and Digital Image Processing, editor 
Kenneth R. Castleman, 1996 Prentice Hall, publisher; The Image Processing 

30 Handbook, editor John C. Russ, 3 rt Edition, 1998, CRC Press; Active 
Contours: The Application of Techniques from Graphics, Vision, Control 
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Theory and Statistics to Visual Tracking of Shapes in Motion, Editors Andrew 
Blake, Michael Isard, 1999 Springer Veriag. As will be appreciated by those of 
skill in the art, as the field of imaging continues to advance methods of 
imaging currently employed can evolve overtime. Thus, any imaging method 
5 or technique that is currently employed is appropriate for application of the 
teachings of this invention as weil as techniques that can be developed in the 
future. A further detailed description of imaging methods is not provided in 
order to avoid obscuring the invention. 

[0050] As shown in Fig. 1a, the first step is to locate a part of the body of a 

10 subject, for example in a human body, for study 98. The part of the body 

located for study is the region of anatomical interest (RAI). In locating a part of 
the body for study, a determination is made to, for example, take an image or 
a series of images of the body at a particular location, e.g. hip, dental, spine, 
etc. Images include, for example, conventional x-ray images, x-ray 

15 tomosynthesis, ultrasound (including A-scan, B-scan and C-scan) computed 
tomography (CT scan), magnetic resonance imaging (MRI), optical coherence 
tomography, single photon emission tomography (SPECT), and positron 
emission tomography, or such other imaging tools that a person of skill in the 
art would find useful in practicing the invention. Once the image is taken, a 

20 region of interest (ROI) can be located within the image 100. Algorithms can 
be used to automatically place regions of interest in a particular image. See, 
e.g., Example 1 describing automatic placement of ROIs in femurs. Image 
data is extracted from the image 102. Finally, quantitative and/or qualitative 
data is extracted from the image data 120. The quantitative and/or qualitative 

25 data extracted from the image includes, for example, the parameters and 
measurements shown in Table 1 , Table 2 or Table 3. 

[0051] Each step of locating a part of the body for study 98, optionally 
locating a region of interest 100, obtaining image data 102, and deriving data 
120, can be repeated one or more times 99,101, 103, 121, respectively, as 
30 desired. 
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[0052] As shown in Fig. 1b image data can be optionally enhanced 104 by 
applying image processing techniques, such as noise filtering or diffusion 
filtering, to facilitate further analysis. Similar to the process shown in Fig. 1a, 
locating a part of the body for study 98, optionally locating a region of 
5 interest 100, obtaining image data 102, enhancing image data 104, and 
deriving data 120, can be repeated one or more times 99,101, 103, 105, 121, 
respectively, as desired. 



Table 1 

Representative Parameters Measured with 
Quantitative and Qualitative Image Analysis Methods 



Parameter 
Bone density and 
nicrostructural 



Measurements """ " ~ " 

'Calibration phantom equivalent thickness 

(Average intensity value of the region of interest expressed as 
thickness of calibration phantom that would produce the equivalent 



■Trabecular contrast 

■Standard deviation of background subtracted ROI 
'Coefficient of Variation of ROI (Standard deviation / mean) 
•(Trabecular equivalent thickness / Marrow equivalent thickness) 
'Fractal dimension 
-Hough transform 
■Fourier spectral analysis 

(Mean transform coefficient absolute value and mean spatial first 



■Predominant orientation of spatial energy spectrum 
Trabecular area 

(Pixel count of extracted trabeculae) 
Trabecular area / Total area 
Trabecular perimeter 

(Count of trabecular pixels with marrow pixels in their neighborhood, 

proximity or vicinity) 
Trabecular distance transform 

(For each trabecular pixel, calculation of distance to closest marrow 

pixel) 

■Marrow distance transform 

(For each marrow pixel, calculation of distance to closest trabecular 
pixel) 

Trabecular distance transform regional maximal values (mean, min., 
max, std. Dev). 

(Describes thickness and thickness variation of trabeculae) 
•Marrow distance transform regional maximal values (mean, min., max, 

std. Dev) 
'Star volume 

(Mean volume of ail the parts of an object which can be seen 
unobscured from a random point inside the object in ail possible 
directions ) 
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'ARAMETER 


Measurements 




•Trabecular Bone Pattern Factor 

(TBPf = (P1 - P2) / (A1 - A2 ) where P1 and A1 are the perimeter 
length and trabecular bone area before dilation and P2 and A2 
corresponding values after a single pixel dilation, measure of 
connectivity) 

•Connected skeleton count or Trees (T) 

•Node count (N) 

•Segment count (S) 

•Node-to-node segment count (NN) 

•Node-to-free-end segment count (NF) 

•Node-to-node segment length (NNL) 

•Node-to-free-end segment length (NFL) 

•Free-end-to-free-end segment length (FFL) 

•Node-to-node total struts length (NN.TSL) 

•Free-end-to-free-ends total struts length( FF.TSL) 

•Total struts length (TSL) 

•FF.TSUTSL 

•NN.TSU TSL 

•Loop count (Lo) 

•Mean distance transform values for each connected skeleton 
•Mean distance transform values for each segment (Tb.Th ) 
•Mean distance transform values for each node-to-node segment 
(Tb.Th. NN) 

•Mean distance transform values for each node-to-free-end segment 

(Tb.Th. NF) 
•Orientation (angle) of each segment 
■Angle between segments 

■Length-thickness ratios (NNL/Tb.Th.NN ) and (NFL/ Tb.Th. NF) 
■Interconnectivity index (ICI) 1CI = (N * NN)/ ( T* (NF + 1) ) 
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Parameter 


Measurements 


Cartilage and cartilage 
tefect/diseased 
:artilage parameters 


•Total cartilage volume 
•Partial/Focal cartilage volume 
•Cartilage thickness distribution (thickness map) 
•Mean cartilage thickness for total region or focal region 
•Median cartilage thickness for total region or focal region 
•Maximum cartilage thickness for total region or focal region 
•Minimum cartilage thickness for total region or focal region 
•3D cartilage surface information for total region or focal region 
•Cartilage curvature analysis for total region or focal region 
•Volume of cartilage defect/diseased cartilage 
•Depth of cartilage defect/diseased cartilage 
•Area of cartilage defect/diseased cartilage 
•2D or '3D location of cartilage defect/diseased cartilage in articular 
surface 

1 «2D or 3D location of cartilage defect/diseased cartilage in 

relationship to weight-bearing area 
•Ratio: diameter of cartilage defect or diseased cartilage / thickness oi 

surrounding normal cartilage 
•Ratio: depth of cartilage defect or diseased cartilage / thickness of 

surrounding normal cartilage 
, «Ratio: volume of cartilage defect or diseased cartilage / thickness of 

surrounding normal cartilage 
•Ratio: surface area of cartilage defect or diseased cartilage / total 

joint or articular surface area 
•Ratio: volume of cartilage defect or diseased cartilage / total cartilage 

volume 


Other articular 
parameters 


•Presence or absence of bone marrow edema 
•Volume of bone marrow edema 

•Volume of bone marrow edema normalized by width, area, size, 
volume of femoral condyle(s)/tibial plateau/patella - other bones 
in other joints 

•Presence or absence of osteophytes 

•Presence or absence of subchondral cysts 

•Presence or absence of subchondral sclerosis 

•Volume of osteophytes 

•Volume of subchondral cysts 

•Volume of subchondral sclerosis 

•Area of bone marrow edema 

•Area of osteophytes 

•Area of subchondral cysts 

•Area of subchondral sclerosis 

•Depth of bone marrow edema 

•Depth of subchondral cysts 
•Depth of subchondral sclerosis 

•Volume, area, depth of osteophytes, subchondral cysts, subchondral 
sclerosis normalized by width, area, size, volume of femoral 
condyte(s)/rjbial plateau/patella -other bones in other joints 

•Presence or absence of menisca! tear 

•Presence or absence of cruciate ligament tear 

•Presence or absence of collateral ligament tear 

•Volume of menisci 

•Ratio of volume of normal to torn/damaged or degenerated meniscal 
tissue 
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Parameter 


Measurements 




•Ratio of surface area of normal to torn/damaged or degenerated 
meniscal tissue 

•Ratio of surface area of normal to torn/damaged or degenerated 
meniscal tissue to total joint or cartilage surface area 

•Ratio of surface area of torn/damaged or degenerated meniscal 
tissue to total joint or cartilage surface area 

•Size ratio of opposing articular surfaces 

•Meniscal subluxation/dislocation in mm 

•Index combining different articular parameters which can also 
include 

oPresence or absence of cruciate or collateral ligament tear 
oBody mass index, weight, height 
•3D surface contour information of subchondral bone 
•Actual or predicted knee flexion angle during gait cycle 
(latter based on gait patterns from subjects with matching 
demographic data retrieved from motion profile database) 
•Predicted knee rotation during gait cycle 
•Predicted knee displacement during gait cycle 
•Predicted load bearing line on cartilage surface during gait cycle and 

measurement of distance between load bearing line and cartilage 

defect/diseased cartilage 
•Predicted load bearing area on cartilage surface during gait cycle 

and measurement of distance between load bearing area and 

cartilage defect/diseased cartilage 
•Predicted load bearing line on cartilage surface during standing or 

different degrees of knee flexion and extension and measurement 

of distance between load bearing line and cartilage 

defect/diseased cartilage 
•Predicted load bearing area on cartilage surface during standing or 

different degrees of knee flexion and extension and measurement 

of distance between load bearing area and cartilage 

defect/diseased cartilage 
•Ratio of load bearing area to area of cartilage defect/diseased 

cartilage 

•Percentage of load bearing area affected by cartilage disease 
•Location of cartilage defect within load bearing area 
•Load applied to cartilage defect, area of diseased cartilage 
•Load applied to cartilage adjacent to cartilage defect area of 
diseased cartilage 



[0053] As will be appreciated by those of skill in the art, the parameters 
and measurements shown in Table 1 are provided for illustration purposes. It 
will be apparent that the terms micro-structural parameters, micro- 

5 architecture, micro-anatomic structure, micro-structural and trabecular 
architecture may be used interchangably. in additon, other parameters and 
measurements, ratios, derived values or indices can be used to extract 
quantitative and/or qualitative information about the ROI without departing 
from the scope of the invention. Additionally, where multiple ROI or multiple 

10 derivatives of data are used, the parameter measured can be the same 
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parameter or a different parameter without departing from the scope of the 
invention. Additionally, data from different ROIs can be combined or 
compared as desired. 

[0054] Additional measurements can be performed that are selected based 
5 on the anatomical structure to be studied as described below. 

[0055] Once the data is extracted from the image it can be manipulated to 
assess the severity of the disease and to determine disease staging (e.g., 
mild, moderate, severe or a numerical value or index). The information can 
also be used to monitor progression of the disease and/or the efficacy of any 
10 interventional steps that have been taken. Finally, the information can be used 
to predict the progression of the disease or to randomize patient groups in 
clinical trials. 

[0056] Fig. 2a illustrates an image 200 taken of an RAI, shown as 202. As 
shown in Fig. 2a, a single region of interest (ROI) 210 has been identified 

15 within the image. The ROI 210 can take up the entire image 200, or nearly the 
entire image. As shown in Fig. 2b more than one ROI can be identified in an 
image. In this example, a first ROI 220 is depicted in one region of the 
image 200, and a second ROI 222 is depicted within the image. In this 
instance, neither of these ROI overlap or abut. As will be appreciated by a 

20 person of skill in the art, the number of ROI identified in an image 200 is not 
limited to the two depicted. Turning now to Fig. 2c another embodiment 
showing two ROI for illustration purposes is shown. In this instance, the first 
ROI 230 and the second ROI 232, are partially overlapping. As will be 
appreciated by those of skill in the art, where multiple ROI are used any or all 

25 of the ROI can be organized such that it does not overlap, it abuts without 
overlapping, it overlaps partially, it overlaps completely (for example where a 
first ROI is located completely within a second identified ROI), and 
combinations thereof. Further the number of ROI per image 200 can range 
from one (ROli) to n (ROI n ) where n is the number of ROI to be analyzed. 

30 [0057] Bone density, microarchitecture, macro-anatomic and/or 

biomechanical (e.g. derived using finite element modeling) analyses can be 
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applied within a region of predefined size and shape and position. This region 
of interest can also be referred to as a "window." Processing can be applied 
repeatedly within the window at different positions of the image. For example, 
a field of sampling points can be generated and the analysis performed at 

5 these points. The results of the analyses for each parameter can be stored in 
a matrix space, e.g., where its position corresponds to the position of the 
sampling point where the analysis occurred, thereby forming a map of the 
spatial distribution of the parameter (a parameter map). The sampling field 
can have regular intervals or irregular intervals with varying density across the 

10 image. The window can have variable size and shape, for example to 
account for different patient size or anatomy. 

[0058] The amount of overlap between the windows can be determined, for 
example, using the interval or density of the sampling points (and resolution of 
the parameter maps). Thus, the density of sampling points is set higher in 

15 regions where higher resolution is desired and set lower where moderate 
resolution is sufficient, in order to improve processing efficiency. The size and 
shape of the window would determine the local specificity of the parameter. 
Window size is preferably set such that it encloses most of the structure being 
measured. Oversized windows are generally avoided to help ensure that 

20 local specificity is not lost. 

[0059] The shape of the window can be varied to have the same 
orientation and/or geometry of the local structure being measured to minimize 
the amount of structure clipping and to maximize local specificity. Thus, both 
2D and/or 3D windows can be used, as well as combinations thereof, 
25 depending on the nature of the image and data to be acquired. 

[0060] In another embodiment, bone density, microarchitecture, macro- 
anatomic and/or biomechanical (e.g. derived using finite element modeling) 
analyses can be applied within a region of predefined size and shape and 
position. The region is generally selected to include most, or all, of the 
30 anatomic region under investigation and, preferably, the parameters can be 
assessed on a pixel-by-pixel basis (e.g., in the case of 2D or 3D images) or a 
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voxel-by-voxel basis in the case of cross-sectional or volumetric images (e.g., 
3D images obtained using MR and/or CT). Alternatively, the analysis can be 
applied to clusters of pixels or voxels wherein the size of the clusters is 
typically selected to represent a compromise between spatial resolution and 
5 processing speed. Each type of analysis can yield a parameter map. 

[0061] Parameter maps can be based on measurement of one or more 
parameters in the image or window; however, parameter maps can also be 
derived using statistical methods. In one embodiment, such statistical 
comparisons can include comparison of data to a reference population, e.g. 
10 using a z-score or a T-score. Thus, parameter maps can include a display of 
z-scores or T-scores. 

[0062] Additional measurements relating to the site to be measured can 
also be taken. For example, measurements can be directed to dental, spine, 
hip, knee or bone cores. Examples of suitable site specific measurements are 
15 shown in Table 2. 
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TABLE 2 

Site specific measurement of bone parameters 



Parameters specific to 
lip images 


•All microarchitecture parameters on structures parallel to stress 
lines 

•All microarchitecture parameters on structures perpendicular to 
stress lines 
Geometry 

•Shaft angle 
•Neck angle 

•Average and minimum diameter of femur neck 
•Hip axis length 

•CCD (caput-collum-diaphysis) angle 

•Width of trochanteric region 

•Largest cross-section of femur head 

•Standard deviation of cortical bone thickness within ROI 

•Minimum, maximum, mean and median thickness of cortical 

bone within ROI 
•Hip joint space width 


Parameters specific to 
spine images 


•All microarchitecture parameters on vertical structures 
•All microarchitecture parameters on horizontal structures 
Seometry 

•Superior endplate cortical thickness (anterior, center, posterior) 
•Inferior endplate cortical thickness (anterior, center, posterior) 
•Anterior vertebral wall cortical thickness (superior, center, 
inferior) 

•Posterior vertebral wall cortical thickness (superior, center, 

•Superior aspect of pedicle cortical thickness 
•inferior aspect of pedicle cortical thickness 
•Vertebral height (anterior, center, posterior) 
•Vertebral diameter (superior, center, inferior), 
•Pedicle thickness (supero-inferior direction). 

•Minimum vertebral height 
•Average vertebral height 
•Anterior vertebral height 
•Medial vertebral height 
•Posterior vertebral height 
•Maximum inter-vertebral height 
•Minimum inter-vertebral height 
•Average inter-vertebral heiqht 


wee images 


•Average medial joint space width 

•Maximum medial joint space width 
•Average lateral joint space width 
•Minimum lateral joint space width 
•Maximum lateral joint space width 



[0063] As will be appreciated by those of skill in the art, measurement and 
5 image processing techniques are adaptable to be applicable to both 
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microarchitecture and macro-anatomical structures. Examples of these 
measurements are shown in Table 3. 



Table 3 

Measurements applicable on Microarchitecture and Macro-anatomical 
5 Structures 



Averaqe density 
measurement 


•Calibrated density of RO! 


Measurements on micro- 
anatomical structures of 
dental, soine, hip, knee or 
aone cores imaoes 


The following parameters are derived from the extracted structures; 
•Calibrated density of extracted structures 
•Calibrated density of background 
•Average intensity of extracted structures 
•Average intensity of background (area other than extracted 
structures) 

•Structural contrast (average intensity of extracted structures / 
average intensity of background ) 

structures / calibrated density of background) 
•Total area of extracted structures 
•Total area of ROI 

•Area of extracted structures normalized by total area of ROI 
•Boundary lengths (perimeter) of extracted normalized by total 
area of ROI 

•Number of structures normalized by area of ROI 
•Trabecular bone pattern factor; measures concavity and 

convexity of structures 
•Star volume of extracted structures 
•Star volume of background 
•Number of loops normalized by area of ROI 


Measurements on 
Distance transform of 
extracted structures 


The following statistics are measured from the distance transform 

eglonal maximum values: 

•Average regional maximum thickness 
•Standard deviation of regional maximum thickness 
•Largest value of regional maximum thickness 
•Median of regional maximum thickness 


skeleton of extracted 
structures 


•Average length of networks (units of connected segments) 
•Maximum length of networks 

•Average thickness of structure units (average distance 

transform values along skeleton) 
•Maximum thickness of structure units (maximum distance 

transform values along skeleton) 
•Number of nodes normalized by ROI area 
•Number of segments normalized by ROI area 
•Number of free-end segments normalized by ROI area 
•Number of inner (node-to-node) segments normalized ROI area 
•Average segment lengths 
•Average free-end segment lengths 
•Average inner segment lengths 
•Average orientation angle of segments 
•Average orientation angle of inner segments 
•Segment tortuosity; a measure of straiqhtness 
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•Segment solidity; another measure of straightness 
•Average thickness of segments (average distance transform 

values along skeleton segments) 
•Average thickness of free-end segments 
•Average thickness of inner segments 
•Ratio of inner segment lengths to inner segment thickness 
•Ratio of free-end segment lengths to free-end segment 

thickness 

•Interconnectivity index; a function of number of inner segments, 
free-end segments and number of networks. 



Directional skeleton 

segment 

measurements 


All measurement of skeleton segments can be constrained by 
one or more desired orientation by measuring only skeleton 
segments within ranges of angle. 


Watershed 
segmentation 


Watershed segmentation is applied to gray level images. 
Statistics of watershed segments are: 
•Total area of segments 

•Number of segments normalized by total area of segments 
•Average area of segments 
•Standard deviation of segment area 
•Smallest segment area 
•Largest segment area 



[0064] As noted above, analysis can also include one or more additional 
techniques include, for example, Hough transform, mean pixel intensity 
analysis, variance of pixel intensity analysis, soft tissue analysis and the like. 
5 See, e.g., co-owned International Application WO 02/30283. 

[0065] Calibrated density typically refers to the measurement of intensity 
values of features in images converted to its actual material density or 
expressed as the density of a reference material whose density is known. The 
reference material can be metal, polymer, plastics, bone, cartilage, etc., and 
10 can be part of the object being imaged or a calibration phantom placed in the 
imaging field of view during image acquisition. 

[0066] Extracted structures typically refer to simplified or amplified 
representations of features derived from images. An example would be binary 
images of trabecular patterns generated by background subtraction and 
1 5 thresholding. Another example would be binary images of cortical bone 

generated by applying an edge filter and thresholding. The binary images can 
be superimposed on gray level images to generate gray level patterns of 
structure of interest. 
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[0067] Distance transform typically refers to an operation applied on binary 
images where maps representing distances of each 0 pixel to the nearest 1 
pixel are generated. Distances can be calculated by the Euclidian magnitude, 
city-block distance, La Place distance or chessboard distance. 

5 [0068] Distance transform of extracted structures typically refer to distance 
transform operation applied to the binary images of extracted structures, such 
as those discussed above with respect to calibrated density. 

[0069] Skeleton of extracted structures typically refer to a binary image of 1 
pixel wide patterns, representing the centerline of extracted structures. It is 
10 generated by applying a skeletonization or medial transform operation, by 
mathematical morphology or other methods, on an image of extracted 
structures. 

[0070] Skeleton segments typically are derived from skeleton of extracted 
structures by performing pixel neighborhood analysis on each skeleton pixel. 
15 This analysis classifies each skeleton pixel as a node pixel or a skeleton 
segment pixel. A node pixel has more than 2 pixels in its 8-neighborhood. A 
skeleton segment is a chain of skeleton segment pixels continuously 8- 
connected. Two skeleton segments are separated by at least one node pixel. 

[0071] Watershed segmentation as it is commonly known to a person of 
20 skill in the art, typically is applied to gray level images to characterize gray 
level continuity of a structure of interest. The statistics of dimensions of 
segments generated by the process are, for example, those listed in Table 3 
above. As will be appreciated by those of skill in the art, however, other 
processes can be used without departing from the scope of the invention. 

25 [0072] Turning now to Fig. 3a, a cross-section of a cartilage defect is shown 
300. The cross-hatched zone 302 corresponds to an area where there is 
cartilage loss. Fig. 3b is a top view of the cartilage defect shown in Fig. 3a. 

[0073] Fig. 3c illustrates the depth of a cartilage defect 310 in a first cross- 
section dimension with a dashed line illustrating a projected location of the 
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original cartilage surface 312. By comparing these two values a ratio of 
cartilage defect depth to cartilage defect width can be calculated. 

[0074] Fig.3d illustrated the depth of the cartilage 320 along with the width 
of the cartilage defect 322. These two values can be compared to determine a 
5 ratio of cartilage depth to cartilage defect width. 

[0075] Fig. 3e shows the depth of the cartilage defect 310 along with the 
depth of the cartilage 320. A dashed line is provided illustrating a projected 
location for the original cartilage surface 312. Similar to the measurements 
made above, ratios between the various measurements can be calculated. 

10 [0076] Turning now to Fig. 3f, an area of bone marrow edema is shown on 
the femur 330 and the tibia 332. The shaded area of edema can be measured 
on a T2-weighted MRI scan. Alternatively, the area can be measured on one 
or more slices. These measurements can then be extended along the entire 
joint using multiple slices or a 3D acquisition. From these measurements 

15 volume can be determined or derived. 

[0077] Fig. 3g shows an area of subchondral sclerosis in the acetabulum 
340 and the femur 342. The sclerosis can be measured on, for example, a T1 
or T2-weighted MRI scan or on a CT scan. The area can be measured on one 
or more slices. Thereafter the measurement can be extended along the entire 
20 joint using multiple slices or a 3D acquisition. From these values a volume can 
be derived of the subchondral sclerosis. For purposes of illustration, a single 
sclerosis has been shown on each surface. However, a person of skill in the 
art will appreciate that more than one sclerosis can occur on a single joint 
surface. 

25 [0078] Fig.3h shows osteophytes on the femur 350 and the tibia 352. The 
osteophytes are shown as cross-hatched areas. Similar to the sclerosis 
shown in Fig. 3g, the osteophytes can be measured on, for example, a T1 or 
T2-weighted MRI scan or on a CT scan. The area can be measured on one or 
more slices. Thereafter the measurement can be extended along the entire 

30 joint using multiple slices or a 3D acquisition. From these values a volume can 
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be derived of the osteophytes. Additionally, a single osteophyte 354 or 
osteophyte groups 356 can be included in any measurement. Persons of skill 
in the art will appreciate that groups can be taken from a single joint surface 
or from opposing joint surfaces, as shown, without departing from the scope 
5 of the invention. 

[0079] Turning now to Fig. 3i an area of subchondral cysts 360, 362, 364 is 
shown. Similar to the sclerosis shown in Fie. 3g, the cysts can be measured 
on, for example, a T1 or T2-weighted MRI scan or on a CT scan. The area 
can be measured on one or more slices. Thereafter the measurement can be 

10 extended along the entire joint using multiple slices or a 3D acquisition. From 
these values a volume can be derived of the cysts. Additionally, single cysts 
366 or groups of cysts 366' can be included in any measurement. Persons of 
skill in the art will appreciate that groups can be taken from a single joint 
surface, as shown, or from opposing joint surfaces without departing from the 

15 scope of the invention. 

[0080] Fig. 3j illustrates an area of torn meniscal tissue (cross-hatched) 
372, 374 as seen from the top 370 and in cross-section 371, Again, similar to 
the sclerosis shown in Fig. 3g, the torn meniscal tissue can be measured on, 
.for example, a T1 or T2-weighted MRI scan or on a CT scan. The area can be 
20 measured on one or more slices. Thereafter the measurement can be 

extended along the entire joint using multiple slices or a 3D acquisition. From 
these values a volume can be derived of the tear. Ratios such as surface or 
volume of torn to normal meniscal tissue can be derived as well as ratios of 
surface of torn meniscus to surface of opposing articulating surface. 

25 [0081] As shown in Fig. 4a, the process of optionally locating a RO1 100, 
extracting image data from the ROI 1 02, and deriving quantitative and/or 
qualitative image data from the extracted image data 120, can be repeated 
122. Alternatively, or in addition, the process of locating a ROI 100, can be 
repeated 124. A person of skill in the art will appreciate that these steps can 

30 be repeated one or more times in any appropriate sequence, as desired, to 
obtain a sufficient amount of quantitative and/or qualitative data on the ROI or 
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to separately extract or evaluate parameters. Further, the ROI used can be 
the same ROI as used in the first process or a newly identified ROI in the 
image. Additionally, as with Fig. 1a the steps of locating a region of interest 
100, obtaining image data 102, and deriving quantitative and/or qualitative 
5 image data can be repeated one or more times, as desired, 101, 103, 121, 
respectively. Although not depicted here, as discussed above with respect to 
Fig. 1a, the additional step of locating a part of the body for study 98 can be 
performed prior to locating a region of interest 100 without departing from the 
invention. Additionally that step can be repeated 99. 

10 [0082] Fig. 4b illustrates the process shown in Fig. 4a with the additional 
step enhancing image data 104. Additionally, the step of enhancing image 
data 104 can be repeated one or more times 705, as desired. The process of 
enhancing image data 104 can be repeated 126 one or more times as 
desired. 

15 [0083] Turning now to Fig. 5a, a process Is shown whereby a region of 
interest is optionally located 100 Although not depicted here, as discussed 
above with respect to Fig. 1a, the step of locating a part of the body for study 
98 can be performed prior to locating a region of interest 100 without 
departing from the invention. Additionally that step can be repeated 99. Once 

20 the region of interest is located 100, and image data is extracted from the ROI 
102, the extracted image data can then be converted to a 2D pattern 130, a 
3D pattern 132 or a 4D pattern 133, for example including velocity or time, to 
facilitate data analyses. Following conversion to 2D 130, 3D 132 or 4D 
pattern 133 the images are evaluated for patterns 140. Additionally images 

25 can be converted from 2D to 3D 131, or from 3D to 4D 131', if desired. 

Although not illustrated to avoid obscuring the figure, persons of skill in the art 
will appreciate that similar conversions can occur between 2D and 4D in this 
process or any process illustrated in this invention. 

[0084] As will be appreciated by those of skill in the art, the conversion 
30 step is optional and the process can proceed directly from extracting image 
data from the ROI 102 to evaluating the data pattern 140 directly 134. 
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Evaluating the data for patterns, includes, for example, performing the 
measurements described in Table 1 .Table 2 or Table 3, above. 

[0085] Additionally, the steps of locating the region of interest 100, 
obtaining image data 102, and evaluating patterns 141 can be performed 
5 once or a plurality of times, 101, 103, 141, respectively at any stage of the 
process. As will be appreciated by those of skill in the art, the steps can be 
repeated. For example, following an evaluation of patterns 140, additional 
image data can be obtained 135, or another region of interest can be located 
137. These steps can be repeated as often as desired, in any combination 
10 desirable to achieve the data analysis desired. 

[0086] Fig. 5b illustrates an alternative process to that shown in Fig. 5a 
which5A that includes the step of enhancing image data 104 prior to 
converting an image or image data to a 2D 130, 3D 132, or 4D 133 pattern. 
The process of enhancing image data 1 04, can be repeated 105 if desired. 
15 Fig. 5c illustrates an alternative embodiment to the process shown in Fig. 5b. 
In this process, the step of enhancing image data 104 occurs after converting 
an image or image data to a 2D 130, 3D 132, or 4D 133 pattern. Again, the 
process of enhancing image data 104, can be repeated 105 if desired. 

[0087] Fig. 5d illustrates an alternative process to that shown in Fig. 5a. 

20 After locating a part of the body for study 98 and imaging, the image is then 
converted to a 2D pattern 130, 3D pattern 132 or 4D pattern 133. The region 
of interest 100 is optionally located within the image after conversion to a 2D, 
3D or 4D image and data is then extracted 102. Patterns are then evaluated 
in the extracted image data 140. As with the process of Fig. 5a, the 

25 conversion step is optional. Further, if desired, images can be converted 
between 2D, 3D 131 and 4D 131' if desired. 

[0088] Similar to Fig. 5a, some or all the processes can be repeated one or 
more times as desired. For example, locating a part of the body for study 98, 
locating a region of interest 100, obtaining image data 102, and evaluating 
30 patterns 140, can be repeated one or more times if desired, 99, 101, 103, 
141, respectively. Again steps can be repeated. For example, following an 
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evaluation of patterns 140, additional image data can be obtained 135, or 
another region of interest can be located 1 37 and/or another portion of the 
body can be located for study 139. These steps can be repeated as often as 
desired, in any combination desirable to achieve the data analysis desired. 

5 [0089] Fig. 5e illustrates an alternative process to that shown in Fig. 5d. In 
this process image data can be enhanced 104. The step of enhancing image 
data can occur prior to conversion 143, prior to locating a region of interest 
145, prior to obtaining image data 102, or prior to evaluating patterns 149. 

[00S0] Similar to Fig. 5a, some or all the processes can be repeated one or 
10 more times as desired, including the process of enhancing image data 104, 
which is shown as 105. 

[0091] The method also comprises obtaining an image of a bone or a joint, 
optionally converting the image to a two-dimensional or three-dimensional or 
four-dimensional pattern, and evaluating the amount or the degree of normal, 

1 5 diseased or abnormal tissue or the degree of degeneration in a region or a 
volume of interest using one or more of the parameters specified in Table 1, 
Table 2 and/or Table 3. By performing this method at an initial time T-i, 
information can be derived that is useful for diagnosing one or more 
conditions or for staging, or determining, the severity of a condition. This 

20 information can also be useful for determining the prognosis of a patient, for 
example with osteoporosis or arthritis. By performing this method at an initial 
time Ti, and a later time T 2 , the change, for example in a region or volume of 
interest, can be determined which then facilitates the evaluation of appropriate 
steps to take for treatment. Moreover, if the subject is already receiving 

25 therapy or if therapy is initiated after time Ti, it is possible to monitor the 
efficacy of treatment. By performing the method at subsequent times, T 2 -Tr,. 
additional data ca be acquired that facilitate predicting the progression of the 
disease as well as the efficacy of any interventional steps that have been 
taken. As will be appreciated by those of skill in the art, subsequent 

30 measurements can be taken at regular time intervals or irregular time 
intervals, or combinations thereof. For example, it can be desirable to 
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perform the analysis at T-i with an initial follow-up, T 2 , measurement taken one 
month later. The pattern of one month follow-up measurements could be 
performed for a year (12 one-month intervals) with subsequent follow-ups 
performed at 6 month intervals and then 12 month intervals. Alternatively, as 
5 an example, three initial measurements could be at one month, followed by a 
single six month follow up which is then followed again by one or more one 
month follow-ups prior to commencing 12 month follow ups. The combinations 
of regular and irregular intervals are endless, and are not discussed further to 
avoid obscuring the invention. 

10 [0092] Moreover, one or more of the parameters listed in Tables 1 , 2 and 3 
can be measured. The measurements can be analyzed separately or the data 
can be combined, for example using statistical methods such as linear 
regression modeling or correlation. Actual and predicted measurements can 
be compared and correlated. See, also, Example 1. 

15 [0093] The method for. assessing the condition of a bone or joint in a 

subject can be fully automated such that the measurements of one or more of 
the parameters specified in Table 1, Table 2 or Table 3 are done 
automatically without intervention. The automatic assessment then can 
include the steps of diagnosis, staging, prognostication or monitoring the 

20 disease or diseases, or to monitor therapy. As will be appreciated by those of 
skill in the art, the fully automated measurement is, for example, possible with 
image processing techniques such as segmentation and registration. This 
process can include, for example, seed growing; thresholding, atlas and 
model based segmentation methods, live wire approaches, active and/or 

25 deformable contour approaches, contour tracking, texture based 

segmentation methods, rigid and non-rigid surface or volume registration, for 
example based on mutual information or other similarity measures. One 
skilled in the art will readily recognize other techniques and methods for fully 
automated assessment of the parameters and measurements specified in 

30 Table 1 , Table 2 and Table 3. 
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[0094] Alternatively, the method of assessing the condition of a bone or 
joint in a subject can be semi-automated such that the measurements of one 
or more of the parameters, such as those specified in Table 1 , are performed 
semi-automatically, i.e., with intervention. The semi-automatic assessment 
5 then allows for human interaction and, for example, quality control, and 
utilizing the measurement of said parameter(s) to diagnose, stage, 
prognosticate or monitor a disease or to monitor a therapy. The semi- 
automated measurement is, for example, possible with image processing 
techniques such as segmentation and registration. This can include seed 

10 growing, thresholding, atlas and model based segmentation methods, live 
wire approaches, active and/or deformable contour approaches, contour 
tracking, texture based segmentation methods, rigid and non-rigid surface or 
volume registration, for example base on mutual information or other similarity 
measures. One skilled in the art will readily recognize other techniques and 

15 methods for semi-automated assessment of the parameters specified in Table 
1, Table 2 or Table 3. 

[0095] Turning now to Fig. 6a, a process is shown whereby the user 
locates a ROI 100, extracts image data from the ROI 102, and then derives 
quantitative and/or qualitative image data from the extracted image data 120, 

20 as shown above with respect to Fig. 1. Following the step of deriving 

quantitative and/or qualitative image data, a candidate agent is administered 
to the patient 150. The candidate agent can be any agent the effects of which 
are to be studied. Agents can include any substance administered or ingested 
by a subject, for example, molecules, pharmaceuticals, biopharmaceuticals, 

25 agropharmaceuticals, or combinations thereof, including cocktails, that are 
thought to affect the quantitative and/or qualitative parameters that can be 
measured in a region of interest. These agents are not limited to those 
intended to treat disease that affects the musculoskeletal system but this 
invention is intended to embrace any and all agents regardless of the 

30 intended treatment site. Thus, appropriate agents are any agents whereby an 
effect can be detected via imaging. The steps of locating a region of interest 
100, obtaining image data 702, obtaining quantitative and/or qualitative data 
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from image data 120, and administering a candidate agent 150, can be 
repeated one or more times as desired, 101, 103, 121, 151, respectively. 

[0096] Fig. 6b shows the additional step of enhancing image data 104, 
which can also be optionally repeated 105 as often as desired. 

5 [0097] As shown in Fig. 6c these steps can be repeated one or more times 
152 to determine the effect of the candidate agent. As will be appreciated by 
those of skill in the art, the step of repeating can occur at the stage of locating 
a region of interest 152 as shown in Fig. 6b or it can occur at the stage 
obtaining image data 153 or obtaining quantitative and/or qualitative data from 
10 image data 154 as shown in Fig. 6d. 

[0098] Fig. 6e shows the additional step of enhancing image data 104, 
which can optionally be repeated 105, as desired. 

[0099] As previously described, some or all the processes shown in 
Figs. 6a-e can be repeated one or more times as desired. For example, 
15 locating a region of interest 100, obtaining image data 102, enhancing image 
data 104, obtaining quantitative and/or qualitative data 120, evaluating 
patterns 140, and administering candidate agent 150 can be repeated one or 
more times if desired, 101, 103, 105, 121, 141, 151 respectively. 

20 [0100] In the scenario described in relation to Figs. 6, an image is taken 
prior to administering the candidate agent. However, as will be appreciated 
by those of skill in the art, it is not always possible to have an image prior to 
administering the candidate agent In those situations, progress is determined 
over time by evaluating the change in parameters from extracted image to 

25 extracted image. 

[0101] Turning now to Fig. 7a, the process is shown whereby the 
candidate agent is administered first 150. Thereafter a region of interest is 
located in an image taken 100 and image data is extracted 102. Once the 
image data is extracted, quantitative and/or qualitative data is extracted from 
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the image data 120. In this scenario, because the candidate agent is 
administered first, the derived quantitative and/or qualitative data derived is 
compared to a database 160 or a subset of the database, which database 
that, includes data for subjects having similar tracked parameters. As shown 
5 in Fig. 7b following the step of obtaining image data, the image data can be 
enhanced 104. This process can optionally be repeated 105, as desired. 

[0102] Alternatively, as shown in Fig. 7c the derived quantitative and/or 
qualitative information can be compared to an image taken at T1 162, or any 
other time, if such image is available. As shown in Fie. 7d the step of 
10 enhancing image data 104 can follow the step of obtaining image data 102. 
Again, the process can be repeated 105, as desired. 

[0103] As previously described, some or all the processes illustrated in 
Figs. 7a-d can be repeated one or more times as desired. For example, 
locating a region of interest 100, obtaining image data 102, enhancing image 

15 data 104, obtaining quantitative and/or qualitative data 120, administering 
candidate agent 150, comparing quantitative and/or qualitative information to 
a database 160, comparing quantitative and/or qualitative information to an 
Image taken at a prior time, such as Ti, 162, monitoring therapy 170, 
monitoring disease progress 172, predicting disease course 174 can be 

20 repeated one or more times if desired, 1 01, 103, 105, 121, 151, 161, 163, 
171, 173, 175 respectively. Each of these steps can be repeated in one or 
more loops as shown in Fig. 7b, 176, 177, 178, 179, 180, as desired or 
appropriate to enhance data collection. 

[01 04] Turning now to Fig. 8a, following the step of extracting image data 
25 from the ROI 102, the image can be transmitted 1 80. Transmission can be to 
another computer in the network or via the World Wide Web to another 
network. Following the step of transmitting the image 180, the image is 
converted to a pattern of normal and diseased tissue 190. Normal tissue 
includes the undamaged tissue located in the body part selected for study. 
30 Diseased tissue includes damaged tissue located in the body part selected for 
study. Diseased tissue can also include, or refer to, a lack of normal tissue in 
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the body part selected for study. For example, damaged or missing cartilage 
would be considered diseased tissue. Once the image is converted, it is 
analyzed 200. Fig. 8b illustrates the process shown in Fig. 8a with the 
additional step of enhancing image data 104. As will be appreciated by those 
5 of skill in the art, this process can be repeated 105 as desired. 

[0105] As shown in Fig. 8c, the step of transmitting the image 180 
illustrated in Fig. 8a is optional and need not be practiced under the invention. 
As will be appreciated by those of skill in the art, the image can also be 
analyzed prior to converting the image to a pattern of normal and diseased. 
10 Fig. 8d illustrates the process shown in Fig. 8c with the additional step of 
enhancing image data 104 that is optionally repeated 105, as desired. 

[0106] As previously described, some or all the processes in Figs. 8a-d 
can be repeated one or more times as desired. For example, locating a region 
of interest 100, obtaining image data 102, enhancing image data 
15 104, transmitting an image 180, converting the image to a pattern of normal 
and diseased 190, analyzing the converted image 200, can be repeated one 
or more times if desired, 101, 103, 105, 181, 191, 201 respectively. 

[0107] Fig. 9 shows two devices 900, 920 that are connected. Either the 
first or second device can develop a degeneration partem from an image of a 

20 region of interest 905. Similarly, either device can house a database for 
generating additional patterns or measurements 915. The first and second 
devices can communicate with each other in the process of analyzing an 
image, developing a degeneration pattern from a region of interest in the 
image, and creating a dataset of patterns or measurements or comparing the 

25 degeneration pattern to a database of patterns or measurements. However, 
all processes can be performed on one or more devices, as desired or 
necessary. 

[0108] In this method the electronically generated, or digitized image or 
portions of the image can be electronically transferred from a transferring 
30 device to a receiving device located distant from the transferring device; 
receiving the transferred image at the distant location; converting the 
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transferred image to a pattern of normal or diseased or abnormal tissue using 
one or more of the parameters specified in Table 1 , Table 2 or Table 3; and 
optionally transmitting the pattern to a site for analysis. As will be appreciated 
by those of skill in the art, the transferring device and receiving device can be 
5 located within the same room or the same building. The devices can be on a 
peer-to-peer network, or an intranet. Alternatively, the devices can be 
separated by large distances and the information can be transferred by any 
suitable means of data transfer, including the World Wide Web and ftp 
protocols. 

10 [0109] Alternatively, the method can comprise electronically transferring an 
electronically-generated image or portions of an image of a bone or a joint 
from a transferring device to a receiving device located distant from the 
transferring device; receiving the transferred image at the distant location; 
converting the transferred image to a degeneration pattern or a pattern of 

15 normal or diseased or abnormal tissue using one or more of the parameters 
specified in Table 1, Table 2 or Table 3; and optionally transmitting the 
degeneration pattern or the pattern of normal or diseased or abnormal tissue 
to a site for analysis. 

[0110] Thus, the invention described herein includes methods and 
20 systems for prognosis of musculoskeletal disease, for example prognosis of 
fracture risk and the like. (See, also, Example 1). Fig. 10 is a schematic 
depiction of an image of a femur showing various ROIs that were analyzed to 
predict fracture risk based on assessment of one or more parameters shown 
in Tables 1 , 2 and 3. 

25 [01 11] In order to make more accurate prognoses, it may be desirable in 
certain instances to compare data obtained from a subject to a reference 
database. For example, when predicting fracture risk, it may be useful to 
compile data of actual (known) fracture load in a variety of samples and store 
the results based on clinical risk factors such as age, sex and weight (or other 

30 characteristics) of the subject from which the sample is obtained. The images 
of these samples are analyzed to obtain parameters shown in Tables 1 , 2 and 
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3. A fracture risk model correlated with fracture load may be developed using 
univariate, bivariate and/or multivariate statistical analysis of these 
parameters and is stored in this database. A fracture risk model may include 
information that is used to estimate fracture risk from parameters shown in 
5 Tables 1, 2 and 3. An example of a fracture risk model is the coefficients of a 
multivariate linear model derived from multivariate linear regression of these 
parameters (Tables 1,2,3, age, sex, weight, etc.) with fracture load. A person 
skilled in the art will appreciate that fracture risk models can be derived using 
other methods such as artificial neural networks and be represented by other 
10 forms such as the coefficients of artificial neural networks. Patient fracture 
risk can then be determined from measurements obtain from bone images by 
referencing to this database. 

[01 12] Methods of determining actual fracture load are known to those in 
the field. Fig. 1 1 is a schematic depiction of biomechanical testing of an intact 
15 femur. As shown, cross-sectional images may be taken throughout testing to 
determine at what load force a fracture occurs. Fig. 12B is a reproduction of 
an x-ray image depicting an example of an induced fracture in a fresh 
cadaveric femur. 

[0113] The analysis techniques described herein can then be applied to a 
20 subject and the risk of fracture (or other disease) predicted using one or more 
of the parameters described herein. As shown in Figs. 13 to 16, the 
prognostication methods described herein are as (or more) accurate than 
known techniques in predicting fracture risk. Fig. 13 is a graph depicting 
linear regression analysis of DXAbone mineral density correlated to fracture 
25 load. Correlations of individual parameters to fracture load are_comparable to 
DXA (Fig. 14 and 15). However, when multiple structural parameters are 
combined, the prediction of load at which fracture will occur is more accurate. 
(Fig. 1 6). Thus, the analyses of images as described herein can be used to 
accurately predict musculoskeletal disease such as fracture risk. 

30 [01 14] Another aspect of the invention is a kit for aiding in assessing the 
condition of a bone or a joint of a subject, which kit comprises a software 
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program, which when installed and executed on a computer reads a 
degeneration pattern or a pattern of normal or diseased or abnormal tissue 
derived using one or more of the parameters specified in Table 1 , Table 2 or 
Table 3 presented in a standard graphics format and produces a computer 
5 readout. The kit can further include a database of measurements for use in 
calibrating or diagnosing the subject. One or more databases can be 
provided to enable the user to compare the results achieved for a specific 
subject against, for example, a wide variety of subjects, or a small subset of 
subjects having characteristics similar to the subject being studied. 

10 [0115] A system is provided that includes (a) a device for electronically 
transferring a degeneration pattern or a pattern of normal, diseased or 
abnormal tissue for the bone or the joint to a receiving device located distant 
from the transferring device; (b) a device for receiving said pattern at the 
remote location; (c) a database accessible at the remote location for 

15 generating additional patterns or measurements for the bone or the joint of the 
human wherein the database includes a collection of subject patterns or data, 
for example of human bones or joints, which patterns or data are organized 
and can be accessed by reference to characteristics such as type of joint, 
' gender, age, height, weight, bone size, type of movement, and distance of 

20 movement; (d) optionally a device for transmitting the correlated pattern back 
to the source of the degeneration pattern or pattern of normal, diseased or 
abnormal tissue. 

[0116] Thus, the methods and systems described herein make use of 
collections of data sets of measurement values, for example measurements of 

25 bone structure and/or bone mineral density from images (e.g., x-ray images). 
Records can be formulated in spreadsheet-like format, for example including 
data attributes such as date of image (x-ray), patient age, sex, weight, current 
medications, geographic location, etc. The database formulations can further 
comprise the calculation of derived or calculated data points from one or more 

30 acquired data points, typically using the parameters listed in Tables 1 , 2 and 
3 or combinations thereof. A variety of derived data points can be useful in 
providing information about individuals or groups during subsequent database 
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manipulation, and are therefore typically included during database 
formulation. Derived data points include, but are not limited to the following: 
(1) maximum value, e.g. bone mineral density, determined for a selected 
region of bone or joint or in multiple samples from the same or different 
5 subjects; (2) minimum value, e.g. bone mineral density, determined for a 
selected region of bone or joint or in multiple samples from the same or 
different subjects; (3) mean value, e.g. bone mineral density, determined for a 
selected region of bone or joint or in multiple samples from the same or 
different subjects; (4) the number of measurements that are abnormally high 

10 or low, determined by comparing a given measurement data point with a 
selected value; and the like. Other derived data points include, but are not 
limited to the following: (1) maximum value of a selected bone structure 
parameter, determined for a selected region of bone or in multiple samples 
from the same or different subjects; (2) minimum value of a selected bone 

15 structure parameter, determined for a selected region of bone or in multiple 
samples from the same or different subjects; (3) mean value of a selected 
bone structure parameter, determined for a selected region of bone or in 
multiple samples from the same or different subjects; (4) the number of bone 
structure measurements that are abnormally high or low, determined by 

20 comparing a given measurement data point with a selected value; and the 
like. Other derived data points will be apparent to persons of ordinary skill in 
the art in light of the teachings of the present specification. The amount of 
available data and data derived from (or arrived at through analysis of) the 
original data provides an unprecedented amount of information that is very 

25 relevant to management of bone-related diseases such as osteoporosis. For 
example, by examining subjects over time, the efficacy of medications can be 
assessed. 

[0117] Measurements and derived data points are collected and 
calculated, respectively, and can be associated with one or more data 
30 attributes to form a database. The amount of available data and data derived 
from (or arrived at through analysis of) the original data provide provides an 
unprecedented amount of information that is very relevant to management of 
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musculoskeletal-related diseases such as osteoporosis or arthritis. For 
example, by examining subjects over time, the efficacy of medications can be 
assessed. 

[0118] Data attributes can be automatically input with the electronic image 
5 and can include, for example, chronological information (e.g., DATE and 
TIME). Other such attributes can include, but are not limited to, the type of 
imager used, scanning information, digitizing information and the like. 
Alternatively, data attributes can be input by the subject and/or operator, for 
example subject identifiers, i.e. characteristics associated with a particular 
10 subject. These identifiers include but are not limited to the following: (1) a 
subject code (e.g., a numeric or alpha-numeric sequence); (2) demographic 
information such as race, gender and age; (3) physical characteristics such as 
weight, height and body mass index (BMI); (4) selected aspects of the 
subject's medical history (e.g., disease states or conditions, etc.); and (5) 
15 disease-associated characteristics such as the type of bone disorder, if any; 
the type of medication used by the subject. In the practice of the present 
invention, each data point would typically be identified with the particular 
subject, as well as the demographic, etc. characteristic of that subject. 

[01 19] Other data attributes will be apparent to persons of ordinary skill in 
20 the art in light of the teachings of the present specification. (See, also, WO 
02/30283, incorporated by reference in its entirety herein). 

[0120] Thus, data (e.g., bone structural information or bone mineral density 
information or articular information) is obtained from normal control subjects 
using the methods described herein. These databases are typically referred 

25 to as "reference databases" and can be used to aid analysis of any given 

subject's image, for example, by comparing the information obtained from the 
subject to the reference database. Generally, the information obtained from 
the normal control subjects will be averaged or otherwise statistically 
manipulated to provide a range of "normal" measurements. Suitable 

30 statistical manipulations and/or evaluations will be apparent to those of skill in 
the art in view of the teachings herein. The comparison of the subject's 
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information to the reference database can be used to determine if the 
subject's bone information falls outside the normal range found in the 
reference database or is statistically significantly different from a normal 
control. 

5 [0121] Data obtained from images, as described above, can be 

manipulated, for example, using a variety of statistical analyses to produce 
useful information. Databases can be created or generated from the data 
collected for an individual, or for a group of individuals, over a defined period 
of time (e.g., days, months or years), from derived data, and from data 
10 attributes. 

[0122] For example, data can be aggregated, sorted, selected, sifted, 
clustered and segregated by means of the attributes associated with the data 
points. A number of data mining software exist which can be used to perform 
the desired manipulations. 

15 [0123] Relationships in various data can be directly queried and/or the data 
analyzed by statistical methods to evaluate the information obtained from 
manipulating the database. 

[01 24] For example, a distribution curve can be established for a selected 
data set, and the mean, median and mode calculated therefor. Further, data 
20 spread characteristics, e.g., variability, quartiles, and standard deviations can 
be calculated. 

[0125] The nature of the relationship between any variables of interest can 
be examined by calculating correlation coefficients. Useful methods for doing 
so include, but are not limited to: Pearson Product Moment Correlation and 
25 Spearman Rank Correlation. Analysis of variance permits testing of 

differences among sample groups to determine whether a selected variable 
has a discernible effect on the parameter being measured. 

[01 26] Non-parametric tests can be used as a means of testing whether 
variations between empirical data and experimental expectancies are 
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attributable to chance or to the variable or variables being examined. These 
include the Chi Square test, the Chi Square Goodness of Fit, the 2x2 
Contingency Table, the Sign Test and the Phi Correlation Coefficient. Other 
tests include z-scores, T-scores or lifetime risk for arthritis, cartilage loss or 
5 osteoporotic fracture. 

[0127] There are numerous tools and analyses available in standard data 
mining software that can be applied to the analyses of the databases that can 
be created according to this invention. Such tools and analysis include, but 
are not limited to, cluster analysis, factor analysis, decision trees, neural 
10 networks, rule induction, data driven modeling, and data visualization. Some 
of the more complex methods of data mining techniques are used to discover 
relationships that are more empirical and data-driven, as opposed to theory 
driven, relationships. 

[0128] Statistical significance can be readily determined by those of skill in 
15 the art. The use of reference databases in the analysis of images facilitates 
that diagnosis, treatment and monitoring of bone conditions such as 
osteoporosis. 

[0129] For a general discussion of statistical methods applied to data 
analysis, see Applied Statistics for Science and Industry, by A. Romano, 
20 1977, Allyn and Bacon, publisher. 

[0130] The data is preferably stored and manipulated using one or more 
computer programs or computer systems. These systems will typically have 
data storage capability {e.g., disk drives, tape storage, optical disks, etc.). 
Further, the computer systems can be networked or can be stand-alone 

25 systems. If networked, the computer system would be able to transfer data to 
any device connected to the networked computer system for example a 
medical doctor or medical care facility using standard e-mail software, a 
central database using database query and update software (e.g., a data 
warehouse of data points, derived data, and data attributes obtained from a 

30 large number of subjects). Alternatively, a user could access from a doctor's 
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office or medical facility, using any computer system with internet access, to 
review historical data that can be useful for determining treatment. 

[0131] If the networked computer system includes a World Wide Web 
application, the application includes the executable code required to generate 
5 database language statements, for example, SQL statements. Such 

executables typically include embedded SQL statements. The application 
further includes a configuration file that contains pointers and addresses to 
the various software entities that are located on the database server in 
addition to the different external and internal databases that are accessed in 
10 response to a user request. The configuration file also directs requests for 
database server resources to the appropriate hardware, as can be necessary 
if the database server is distributed over two or more different computers. 

[0132] As a person of skill in the art will appreciate, one or more of the 
parameters specified in Table 1, Table and Table 3 can be used at an initial - 
15 time point Ti to assess the severity of a bone disease such as osteoporosis or 
arthritis. The patient can then serve as their own control at a later time point 
T 2 , when a subsequent measurement using one or more of the same 
parameters used at Ti is repeated. 

[0133] A variety of data comparisons can be made that will facilitate drug 
20 discovery, efficacy, dosing, and comparisons. For example, one or more of 
the parameters specified in Table 1, Table 2 and Table 3 may be used to 
identify lead compounds during drug discovery. For example, different 
compounds can be tested in animal studies and the lead compounds with 
regard to highest therapeutic efficacy and lowest toxicity, e.g. to the bone or 
25 the cartilage, can be identified. Similar studies can be performed in human 
subjects, e.g. FDA phase I, II or III trials. Alternatively, or in addition, one or 
more of the parameters specified in Table 1, Table 2 and Table 3 can be used 
to establish optimal dosing of a new compound. It will be appreciated also that 
one or more of the parameters specified in Table 1, Table 2 and Table 3 can 
30 be used to compare a new drug against one or more established drugs or a 
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placebo. The patient can then serve as their own control at a later time point 
T 2l 

EXAMPLES 

Example 1: Correlation of Macro-anatomical and Structural 
5 Parameters to Fracture Load 

[0134] Using 15 fresh cadaveric femurs, the following analyses were 
performed to determine the correlation of macro-anatomical and structural 
parameters to fracture load. 

[0135] Standardization of Hip radiographs : Density and magnification 
10 calibration on the x-ray radiographs was achieved using a calibration 
phantom. The reference orientation of the hip x-rays was the average 
orientation of the femoral shaft. 

[0136] Automatic Placement of Regions of Interest. An algorithm was 
developed and used to consistently and accurately place 7 regions of interest 

15 based on the geometric and position of proximal femur. Fig. 10. In brief, the 
algorithm involved the detection of femoral boundaries, estimation of shaft 
and neck axes, and construction of ROI based on axes and boundary 
intercept points. This approach ensured that the size and shape of ROIs 
placed conformed to the scale and shape of the femur, and thus were 

20 consistent relative to anatomic features on the femur. 

[0137] Automatic Segmentation of the proximal femur. A global gray level 
thresholding using bi-modal histogram segmentation algorithm(s) was 
performed on the hip images and a binary image of the proximal femur was 
generated. Edge-detection analysis was also performed on the hip x-rays, 
25 including edge detection of the outline of the proximal femur that involved 
breaking edges detected into segments and characterizing the orientation of 
each segment. Each edge segment was then referenced to a map of 
expected proximal femur edge orientation and to a map of the probability of 
edge location. Edge segments that did not conform to the expected 
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orientation or which were in low probability regions were removed. 
Morphology operations were applied to the edge image(s) to connect any 
discontinuities. The edge image formed an enclosed boundary of the 
proximal femur. The region within the boundary was then combined with the 
5 binary image from global thresholding to form the final mask of the proximal 
femur. 

[0138] Automatic Segmentation and Measurement of the Femoral Cortex : 
Within a region of interest (ROI), edge detection was applied. Morphology 
operations were applied to connect edge discontinuities. Segments were 

10 formed within enclosed edges. The area and the major axis length of each 
segment were then measured. The regions were also superimposed on the 
original gray level image and average gray level within each region was 
measured. The cortex was identified as those segments connected to the 
boundary of the proximal femur mask with the greatest area, longest major 

15 axis length and a mean gray level about the average gray level of ail enclosed 
segments within the proximal femur mask. 

[0139] The segment identified as cortex was then skeletonized. The 
orientation of the cortex skeleton was verified to conform to the orientation 
map of the proximal femur edge. Euclidean distance transform was applied to 
20 the binary image of the segment. The values of distance transform value 
along the skeleton were sampled and their average, standard deviation, 
minimum, maximum and mod determined. 

[0140] Watershed Segmentation for Characterizing Trabecular Structure : 
Marrow spacing was characterized by determining watershed segmentation of 

25 gray level trabecular structures on the hip images, essentially as described in 
Russ "The Image Processing Handbook," 3 rd . ed. pp.494-501. This analysis 
take the gray level contrast between the marrow spacing and adjacent 
trabecular structures into account. The segments of marrow spacing 
generated using watershed segmentation were measured for the area, 

30 eccentricity, orientation, and the average gray level on the x-ray image within 
the segment. Mean, standard deviation, minimum, maximum and mod. were 
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determined for each segment. In addition, various structural and/or macro- 
anatomical parameters were assessed for several ROIs (FIG. 10). 

[0141] Measurement of Femoral Neck BMP;. DXA analysis of bone mineral 
density was performed in the femoral neck region of the femurs. 

5 [01421 Biomechanical Testing of Intact Femur . Each cadaveric femur 
sample (n=1 5) was tested for fracture load as follows. First, the femur was 
placed at a 15° angle of tilt and an 8° external rotation in an Instron 1331 
Instrument (Instron, Inc.) and a load vector at the femoral head simulating 
single-leg stance was generated, essentially as described in Cheal et al. 

10 (1992) J. Orthop. Res. 10(3):405-422. Second, varus/ valgus and torsional 
resistive movements simulating passive knee ligaments restraints were 
applied. Next, forces and movement at failure were measured using a six- 
degree of freedom load cell. Subsequently, a single ramp, axial compressive 
load was applied to the femoral head of each sample at 100 mm/s until 

15 fracture. (FIG. 12). Fracture load and resultant equilibrium forces and 

moments at the distal end of the femur were measured continuously. FIG. 1 1 
shows various results of biomechanical testing. 

[0143] The correlation between (1) DXA femoral next BMD and facture 
load; (2) bone structure and fracture load; and (3) macro-anatomical analyses 
20 and fracture load was determined and shown in FIG. 13-15, respectively. 

[0144] Multivariate linear regression analysis was also performed, 
combining several structural and macro-anatomical parameters, including 
local maximum marrow spacing (r=0.6 linearized); standard deviation of 
cortical thickness of ROI3 (r=0.57); maximum cortical thickness of ROI5 

25 (r=o!56); and mean node-free end length for ROI3 (r=0.50). Results are 
shown in FIG. 16 and demonstrate that, using analyses, described herein 
there is a good correlation between predicted fracture load and actual fracture 
load (r=0.81, p<0.001). The mean fracture load was 5.4 kiloNewton with a 
standard deviation of 2.3 kiloNewton. These statistics and the coefficients of 

30 multivariate linear regression were stored as data of the fracture load 
reference database. 
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Example 2: Correlation of 2D and 3D measurements 

[0145] To demonstrate that methods using 2D x-ray technology to 
quantitatively assess trabecular architecture is as effective as 3D nCT, which 
serves as a gold standard for such measurements, the following experiments 
5 were performed. Bone cores (n=48) were harvested from cadaveric proximal 
femora. Specimen radiographs were obtained and 2D structural parameters 
were measured on the radiographs. Cores were then subjected to 3D nCT 
and biomechanical testing. The nCT images were analyzed to obtained 3D 
micro-structural measurements. Digitized 2D x-ray images of these cores 
10 were also analyzed as described herein to obtain comparative micro-structural 
measurements. 

[0146] Results showed very good correlation among the numerous 2D 
parameters and 3D (iCT measurements, including for example correlation 
between 2D Trabecular Perimeter/Trabecular Area (Tb.P/Tb.A) with 3D Bone 

15 Surface/Bone Volume (r = 0.92, p< 0.001), and 2D Trabecular Separation 
(Tb.Sp) with 3D Trabecular Separation (r = 0.88, p < 0.001). The 2D 
Tb.P/Tb.A and 2D Tb.Sp also function correlate very well as predictive 
parameters for the mechanical loads required to fracture the cores, with r = - 
0.84 (p <0.001) and r = -0.83 (p < 0.001), respectively, when logarithmic and 

20 exponential transformations were used in the regression, 

[0147] These results demonstrate that 2D micro-structural measurements of 
trabecular bone from digitized radiographs are highly correlated with 3D 
measurements obtained from n-CT images. Therefore, the mechanical 
characteristics of trabecular bone microstructure from digitized radiographic 
25 images can be accurately determined from 2D images. 

Example 3. Prediction of fracture risk using fracture load 
reference database. 

[0148] A hip x-ray of cadaver pelvis was exposed using standard clinical 
procedure and equipment. The radiograph film was developed and digitized. 
30 The image was then analyzed to obtain micro-structure, and macro- 
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anatomical parameters. The local maximum spacing, standard deviation of 
cortical thickness of ROI3, maximum cortical thickness of ROI5, and mean 
node-free end length for ROI3 were used to predict load required to fracture 
the cadaver hip using the coefficients of multivariate linear regression stored 
5 in the fracture load reference database. The predicted fracture load was 7.5 
kiloNewton. This fracture load is 0.98 standard deviation above the average of 
the fracture load reference database (or z-score = 0.98). This result may 
suggest that the subject had a relatively low risk of sustaining a hip fracture as 
compared to the population of the reference database. 

10 [0149] The foregoing description of embodiments of the present invention 
has been provided for the purposes of illustration and description. It is not 
intended to be exhaustive or to limit the invention to the precise forms 
disclosed. Many modifications and variations will be apparent to the 
practitioner skilled in the art. The embodiments were chosen and described in 

15 order to best explain the principles of the invention and its practical 
application, thereby enabling others skilled in the art to understand the 
invention and the various embodiments and with various modifications that 
are suited to the particular use contemplated. It is intended that the scope of 
the invention be defined by the following claims and its equivalence. 



47 



WO 2004/062495 



PCT/US2004/000315 



What is claimed is: 

1. A method of predicting bone or articular disease in a subject, the 
method comprising the steps of: 

5 determining one or more micro-structural parameters, one or more 

macroanatomical parameters or biomechanical parameters of a joint in said 
subject; and 

combining at least two of said parameters to predict the risk of bone or 
articular disease. 

10 

2. The method of claim 1 , wherein said combining comprises 
combining one or more micro-structural parameters and one or more macro- 
anatomical parameters. 

15 3. The method of claim 1, wherein said combining comprises 

combining one or more micro-structural parameters and one or more 
biomechanical parameters. 

4. The method of claim 1, wherein said combining comprises 
20 combining one or more macroanatomical parameters and one or more 

biomechanical parameters. 

5. The method of claim 1 , wherein said combining comprises 
combining one or more macroanatomical parameters, one or more micro- 

25 structural parameters and one or more biomechanical parameters. 

6. The method of any of claims 1 to 5, wherein said bone or articular 
disease is fracture risk. 

30 7. The method of any of claims 1 to 5, wherein the parameters are 

obtained from one or more regions of interest in an image obtained from said 
subject. 
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8. The method of claim 7, wherein the image comprises a calibration 
phantom. 

9. The method of any of claims 1 to 5, wherein said parameters are 

5 selected from the group consisting of one or more of the parameters set forth 
in Tables 1, 2 and 3. 

10. The method of any of claims 1 to 5, wherein said combining 
comprises univariate, bivariate or multivariate statistical analysis. [ 

10 

11. The method of any of claims 1 to 5, further comprising comparing 
said parameters to data derived from a reference database of known disease 
parameters. 

15 12. The method of any of claims 1 to 5, wherein the bone is in a region 

selected from the group consisting of leg, knee, hip, spine and arm. 

13. The method of claim 7, wherein the image is selected from the 
group consisting of an x-ray image, a CT image, an ultrasound image and an 

20 MRI. 

14. The method of any of claims 1 to 5, further comprising 
administering a compound to the subject. 

25 15. The method of claim 14, wherein the steps are repeated at two or 

more time points and further wherein one time point is prior to administration 
of the compound. 

16. A method of determining the effect of a candidate agent on a 
30 subject's prognosis for musculoskeletal disease comprising: 

predicting a first risk of musculoskeletal disease in subject according to 
the method of any of claims 1 to 5; 

administering a candidate agent to said subject; 
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predicting a second risk of said musculoskeletal disease in said subject 
according to the method of any of claims 1 to 5; and 

comparing said first and second risks, thereby determining the effect of 
the candidate on the subject's prognosis for said disease. 

5 ■ 

17. The method of claim 16, wherein said candidate agent is 
administered to the subject. 

18. The method of claim 16, wherein said administration comprises 
10 ingestion or injection. 

19. The method of claim 16, wherein said candidate agent is selected 
from the group consisting of molecules, pharmaceuticals, biopharmaceuticals, 
agropharmaceuticals and combinations thereof. 

15 
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f. Q schemes, rules or methods of doing business. 

g. Q schemes, rules or methods of performing purely mental acts. 

h. Q schemes, rules or methods of playing games. 

i. \x\ nsthodsfortreatmentofthehumanbodybysurgeryortherapy. 
j. [x] methods for treatment of the animal body by surgery or therapy, 
k- [x] diagnostic methods practised on the human or animal body. 

I. Q mere presentations of information. '"' 

m. Q computer programs lor which this International Searching Authority is not equipped to 



□ the description [x]. the claims Q the drawings 

3. □ The failure of the nucleotide and/or amino acid sequence listing to comply with the standard provided for in Annex C of the 

Administrative Instructions prevents a meaningful search from being carried out: 

f~1 the written form has not been furnished or does not complywith the standard. 

Q the computer readable form has not been furnished or does not comply with the standard: 

4. Q The failure of the tables related to the nucleotide and/or amino acid sequence listing to comply with the technical 

- requirements provided for In Annex C-Ws of the Administrative Instructions prevents a meaningful search from being 
carried out: 

I | the written form has not been furnished. 

□ the computer readable form has not been furnished or does not comply with the technical requirements. 
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FURTHER INFORMATION CONTINUED FROM PCT/ISA/ 203 



A meaningful search is not possible on the basis of all claims because 
all claims are directed to - Diagnostic method practised on the human or 
animal body - Rule 39.1(iv) PCT for claims 1-19 

Furthermore, claims 14-17 refer to method for treatment of the human body 
by therapy - Rule 39.1(iv) 

And claim 18 also refers to method for treatment of the human body' by 
surgery - Rule 39.1 (iv) 

The applicant's attention is drawn to the fact that claims relating to 
inventions in respect of which no international search report has been 
established need not be the subject of an international preliminary 
examination (Rule 66.1(e) PCT). The applicant is advised that the EPO 
policy when acting as an International Preliminary Examining Authority is 
normally not to carry out a preliminary examination on matter which has 
not been searched. This is the case irrespective of whether or not the 
claims are amended following receipt of the search report or during- any 
Chapter II procedure. If the application proceeds into the regional phase 
before the EPO, the applicant is reminded that a search may be carried 
out during examination before the EPO (see EPO Guideline C-VI, 8.5), 
should the problems which led to the Article 17(2) declaration be 
overcome. 



